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Abstract 



We describe a measurement of B°—B° mixing parameters exploiting a method 
of partial reconstruction of the decay chains B — > D* ± 7r =F and B — > D*^p T . 
Using 9.6 xlO 6 BB pairs collected at the Cornell Electron Storage Ring, we 
find Xd = 0.198 ± 0.013 ± 0.014, \y d \ < 0.41 at 95% confidence level, and 
\Ue(e B )\ < 0.034 at 95% confidence level. 
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The discovery of B° — B° mixing in 1987 [l] signaled a large top quark mass and allows 
the anticipation of observable CP violating asymmetries in the B° meson in the near future. 
Well-known values of the parameters describing mixing will be necessary to extract precise 
values of CP violating parameters, as planned at asymmetric B factories 0. In addition, 
the size of B° — B° mixing is characterized by the mass difference parameter Arrid and is 
proportional to the square of |V^V^|, the magnitude of one side of the "Unitarity Triangle" 
which describes some of the mathematical constraints imposed by unitarity upon the 
elements of the CKM matrix Accurate measurements of Am for B° and B° s mesons 
therefore provide an independent check on our understanding of CP violation in the Standard 
Model. 

Mixing in the B° — B° system may be described by the parameters Xd, yd, P, and q. 
The parameter Xd = Arrid /Td, where Arrid is the mass difference between the heavy and 
light eigenstates B H and B L and T d is the average natural decay width. Similarly, y d is the 
normalized lifetime difference between the two eigenstates, and can be written as ATd/2Td- 
The parameters p and q describe the B° and B° amplitudes, respectively, in the eigenstates 
Bh and B^. When the T(4S) is produced with symmetric energy electron positron collisions, 
the experimentally accessible quantity is Xd, the time-integrated probability for an initially 
produced B° or B° meson to decay as its CP conjugate. It may be written in terms of Xd 
and y d : 

= T(B°^W) _ xl + yj 

Xd T(B° -> B°) + T(B° -f B°) 2(1 + 3%)' {) 



Under certain assumptions, a measurement of Xd can be combined with direct determi- 
nations of Arrid and the B meson lifetime in order to extract yd- The ancillary vari- 
able, 6b, is analogous to the K°- mixing parameter e, and is defined through the relation 
p = (1 + e#) / \J2{1 + |ee| 2 ). Limits on e B can be extracted by searching for a CP violating 
asymmetry in the events where mixing has occurred. The mixing parameters yd and e# are 
both expected to be of order 10~ 2 with considerable uncertainty. 

In this letter, we report new measurements of the B° — B° mixing parameters measured 
at the T(4S) resonance. We attempt to determine the beauty quantum number (or flavor) at 
decay for both of the B mesons produced in the T(4S) decay using a novel method subject 
to systematic uncertainties very different from previous measurements 0. When the 
decay flavors of the two B mesons in the event coincide, it indicates that the second B has 
undergone mixing in the interval between the decays of the two B mesons. The flavor of 
one B meson at decay is tagged by a high-momentum lepton originating from the decay 
chain B — > Xlv. The flavor at decay of the remaining B meson is determined through 
partial reconstruction of the decay chain B° — > D* + h^, (charge conjugate modes implied), 
where refers either to a 7r~ or p~. The electric charge of the h\y identifies (tags) the 
value of the B flavor at the time of its decay. (We assume the double Cabibbo suppressed 
decay B° — > D*~h^ is negligible.) By employing the hadronic flavor tag for one B° in the 
event, this method sacrifices statistical accuracy relative to methods where a semileptonic 
decay is used to tag the flavor of both B mesons in the event. However, the systematic error 
due to the uncertainty in the charged to neutral B meson production ratio at the T(4S") 
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that dominates measurements of Xd at the T(4S) using dileptons is substantially reduced 
||. As a result this method results in a significant improvement in precision over previous 
measurements of Xd at the T(4S f ). 

Four charge combinations of hadrons and leptons are possible: h w £ + , h w £~ , h w £~ , 
and h w £ + . In the absence of backgrounds or mistags, these correspond to the four flavor 
combinations B°B°, WW, B°W, and WB°, respectively. Then, 

h w £ + + h w £ . . 

Xd h+£ + + hZri- + hU~ + ' 1 1 



In practice, the raw counts recorded in each of the combinations must be corrected for 
processes that incorrectly tag the B decay flavor (mistags). Mistags may be due either 
to leptons not arising from the primary decay B — > X£v or hadrons not arising from the 
hypothesized decay chain B° — > D* + h^ v . Backgrounds (non-i? events), which can contribute 
either to the denominator or numerator, must be subtracted. 

The data were recorded at the Cornell Electron Storage Ring (CESR) with two con- 
figurations of the CLEO detector called CLEO II and CLEO II. V. In the CLEO II. V 
configuration, the innermost wire chamber was replaced with a precision three-layer silicon 
vertex detector (SVX) ||. The results presented here are based upon an integrated lumi- 
nosity of 9.1 fb _1 of e + e~ data taken at the T(45) energy and 4.4 fb _1 taken an average of 
60 MeV below BB threshold. The Monte Carlo simulation of the CLEO detector was based 
upon GEANT || and simulated events were processed in the same manner as the data. 

The method of partial reconstruction used in this Letter has been described in detail 
elsewhere ||10|| . By observing the hw and the soft pion from the decay D* + — ► D tt^, we 
deduce the kinematics of the decay chain B° — > D* + h w without reconstruction of the D°. 
We denote the charged pion from the h w as itf. 

The hadronic B° decay may be described by three angles. The angle formed, in the B 
(D*) rest frame, between the D*(D°) flight direction and the direction of the lab frame, 
is called 0* B {6* D *). A larger value of cos 6* B (cos 9* D *) corresponds to a higher momentum of 
the hw(n s ). The third angle, 0, is the angle between the plane of the B° — > D* + h w decay 
and the plane that contains the h w and the 7r+, as shown in Fig. ffl. All three angles have 
distinctive distributions for signal and background. The cos9 B distribution is constant for 
signal because the B meson is a scalar particle. The distribution of cos 6^* shows the 100% 
polarization in the D* + tt~ mode from conservation of angular momentum, and shows the 
87% polarization that has been measured in the D* + p~ mode. The distribution of cos<ft is a 
combination of the ir s and hw momenta and the angle between them. For most signal events 
it reconstructs inside the physical region — 1 < cos0 < 1. For non-signal events as well as 
signal events with imperfect measurements of the pion momenta, cos <fi can be calculated but 
may fail to have a physical value. 

In order to select hadronic D* + h^ decay candidates for the analysis, both the ir s and 
iTf candidate tracks must be well-reconstructed and consistent with originating at the e + e~ 
interaction point, and must not be identified as a lepton. We reconstruct candidates from 
7Tj 7T° combinations, where the 7r° is formed from a pair of photon candidates. 

The hw momentum is required to fall in the kinematically allowed range for B° — > D* + h w 
decays, assuming E B o = Et, eam . We require the momentum of the n s to be below 300 MeV/c. 
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The tt s and 7Tj are required to have opposite electrical charges. We require | cos0| < 7 and 
cos 0h w -K S < —0.8, where 9h w n a is the angle between the h w and n s . In 8% of selected events, 
there is more than one combination of charged tracks that satisfy these criteria. In this case, 
we select the one that is reconstructed as D*it rather than D*p. If more than one combination 
still remains, we choose the one for which the value of cos0 is nearest 0.6 (the peak of the 
signal distribution.) The resulting bias in the cos</> shape has a negligible effect upon our 
measurement of the mixing parameters. The events satisfying these criteria are used to 
determine the sample composition. For the events we use to measure yields (and therefore the 
mixing parameters), we use more restrictive criteria, requiring cos hw7Ts < —0.95, | cos <f>\ < 2, 



and | cos9* D , | < 1. 

Lepton candidates are selected by requiring that the track is well-reconstructed, consis- 
tent with originating from the e + e~ interaction point, and well identified as either an electron 
or a muon. We require that the momentum of the lepton candidate is greater than 1.4 GeV/c 
and we use the angle between the lepton and the n s to suppress semileptonic decays from 
the unreconstructed D°. We veto leptons from the decays J/tp — > £ + £~. If the event has 
been partially reconstructed as D* + n~, we require that the lepton form a large angle with 
the thrust axis of the remainder of the event in order to suppress e + e~ — > qq backgrounds, 
where q = {u,d, s,c}. If more than one lepton candidate in an event satisfies the criteria, 
we select the highest momentum candidate. 

For the events selected by these criteria, the B decay modes contributing to the D* + h : w 
candidates can be divided into five categories: signal (D* + n~ and D* + p~), other two-body 
and semi-leptonic B° decays (such as B° — > D + p~ and B° — > D + £~V), two-body and 
semi-leptonic B^ decays (such as B~ — > D*°tt~ and B~ — > D*°£~V), random combinatoric 
backgrounds, and events of the type e + e~ — > qq (continuum). For the distributions of 
two-body B° and B + decays in cos9* D * and cos</>, we rely on the simulation. We include 
10 two-body and semi-leptonic decay modes of the B° in the definition of two-body B° 
decays, and 12 in the definition of the two-body B + decays |J2]. These decays are well- 
measured H and are reliably modeled by the simulation. Combinations of h\y and tc s that 
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satisfy the analysis requirements, yet originate from neither the signal decays nor the two- 
body decays, are considered random combinatoric backgrounds. To model these, we use 
a synthetic distribution in cosflj,. and cos0, generated by combining track pairs drawn 
at random from the observed spectrum of all BB track momenta in data. The cosine of 
the angle between them, cos 6^^, is distributed uniformly. The simulation predicts that 
the distribution generated with this procedure provides an excellent approximation to the 
distributions of these decays, and indicates that 40% of the combinatoric background comes 
from B + B~ events. Distributions from continuum qq production are directly measured in 
data taken below the T(AS) resonance. 

In order to determine the composition of our event sample, we divide the sample into 
two subsets based on the value of cos 9h wns - Events for which —0.9 < cos 9h wlTa < —0.8 are 
the sideband sample; events for which cos 9h w7T3 < —0.95 make up the signal sample. We 
perform a binned two-dimensional maximum likelihood fit jO] simultaneously to the cos(/> 
and cos 9* Di( distributions of the signal sample and the sideband sample. The fit determines 
the normalization of the first four categories of events; that of the continuum is fixed by the 
relative luminosity of the data taken at and below the T(4S). 

The results of the fit are shown in Table |. The two projections of the fits in the mixing 
sample (more restrictive selection criteria) are shown in Fig. ^|. We find that the fit has a 
Baker-Cousins \ 2 IHI of 129.7 for 151 degrees of freedom. From simulations we find that 



this corresponds to a confidence level of 91%. 
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FIG. 2. Projections of the fit to the signal sample onto (a) the cos<^> axis and (b) the cos8d* 
axis. In each plot, the points are the data and the histograms are the best-fit shapes from the 
simulation. 



We also show in Fig. that the distributions of the subset of events in Fig. ^| that contribute 
to the numerator of Eqn. [| (h^l^ 1 events), are well described by the fit. 

The partially reconstructed hadronic tag may incorrectly identify the flavor of the de- 
caying B meson. The dominant source of mistagged events is B° candidates formed from 
random combinations of tracks in B + B~ or B°B° events. We determine that the mistag 
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FIG. 3. Projections of the fit to the sample of mixed events onto (a) the cos<^> axis and (b) the 
cos#d* axis. In each plot, the points are the data and the histograms are the best-fit shapes from 
the simulation. 

TABLE I. The composition of the mixing sample, as determined by the fit. For comparison, 
the continuum subtracted data sample size is 1865 events. The mistag fraction uncertainties are 
statistical only. 



Event Type 



Fit 



Mistag Fraction 



Signal 
Two-body B° 
Two-body 
Combinatoric 

Total 



1241±52 
262±60 
172±45 
192±21 
1867±92 



0.0006±0.0006 
0.020±0.005 
0.050±0.010 

0.21±0.12 
0.031±0.012 



rate of combinatoric events is (21±12)% (where the uncertainty is statistical only), using 
a separate sample of fully-reconstructed B — > D*tv decays in the data. We combine the 
composition as determined by the fit (Table |) with the mistag rates for each individual 
component to determine the hadronic mistag rate, also shown in Table |I[ We calculate a to- 
tal hadronic mistag rate of (3.1±1.2)%. The uncertainty includes the statistical uncertainties 
from the fit and in the random combinatoric mistag rate. 

The lepton in the event may also mistag the flavor of the decaying B meson for several 
reasons. Leptons may arise from the secondary decay chain B — > DX, D — > Xtv. The 
magnitude of this source is well-constrained by measurements of the B — > DX spectrum 
]14| and the known form factors governing D semileptonic decays []15 |. We also correct for 



leptons from qq events, B — > ipX, ip ~* ^ + ^~ events ||, misidentified hadrons fll2"| , leptons 
from [IB] and other upper vertex (6 — > c) production |T^], in-flight decays, it Dalitz 
decays, 7 conversions, and 5 rays [JjTJ]. Altogether we find that (3.6±0.5)% of electrons and 
(3.8±0.5)% of muons incorrectly tag the B decay flavor. The uncertainties are the total 
systematic uncertainties obtained by adding in quadrature the uncertainties associated with 
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TABLE II. Data event yields, sequentially subjected to the corrections, and the value of Xd 
computer from these yields. The uncertainties are statistical only. 



Correction 



h 



±0± 



h ± £* 



Corrected x 



None (raw yield) 

Continuum 
Electron mistags 
Muon mistags 
Hadron mistags 
B + B~ Background 



458±21 
401±23 
377±23 
359±24 
321±25 
321±25 



1524±39 
1464±40 
1487±31 
1505±42 
1543±43 
1294±38 



0.231±0.010 
0.215±0.011 
0.202±0.011 
0.193±0.012 
0.172±0.012 
0.198±0.013 



TABLE III. Summary of systematic uncertainties for Xd- 



Source Uncertainty 

Hadronic Mistag Fraction 0.009 

B ± Background 0.007 

Two-Body Distributions 0.006 

Mixed:Unmixed Efficiency Difference 0.004 

Lepton Mistag Fraction 0.003 

Mixed: Unmixed Mistag Rate Difference 0.003 

Signal Shape 0.002 

Combinatoric Shape 0.002 

Total 0.014 



each of the input branching fractions, spectra, and fake rates. 

The yields for the mixing sample (more restrictive selection criteria) in the possible 
charge combinations, and subsequent corrections, are summarized in Table We correct 
the continuum subtracted raw yields of Table [TT] for the mistag levels that we have determined 
for the leptonic and hadronic tags, then subtract B + B~ background, which contributes to 
the denominator even when the beauty quantum number is correctly reconstructed. The 
total charged B background is (13.3±2.5)%, which is the sum of the two-body B + B~ decays 
and the 40% of the combinatoric background that is attributed to B + B~ . The fully-corrected 



result is Xd = 0.198 ± 0.013 ± 0.014. The systematic uncertainties are listed in Table [II 



The largest systematic uncertainty in Xd is due to the uncertainty in the total hadronic 
mistag rate which in turn is dominated by the uncertainty in the mistag rate of combinatoric 
background events. The systematic uncertainty in the combinatoric background mistag rate 
is determined by comparing the mistag rates of energetic pions in data and simulated events, 
using samples in which the B flavor has been tagged using the decay D* + £~T7. Smaller 
contributions to the total mistag rate come from the statistical uncertainty of the fit and 
uncertainty in the two-body mistag rates. We evaluate mistag rates for hadronic and leptonic 
tags independently, assuming no correlation. The difference in efficiencies and mistag rates 
for mixed and unmixed events are both found to be consistent with zero in large samples of 
simulated signal events, and the uncertainty in Xd reflects the statistical uncertainty of the 
finite simulation samples. We assign a systematic uncertainty in Xd due to the uncertainty 
in the mistag rate totaling 9%. 
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The uncertainty in the B + B~~ background is dominated by uncertainty in the percentage 
of random combinations that are due to B + B~ decays. The assigned uncertainty allows 
the fraction of random combinations arising from B + B~ decays to vary uniformly from to 
100%. 

The uncertainty due to the distributions in cos 9* D * and cos <fi of two-body decays is evalu- 
ated by repeating the analysis, varying in turn each two-body decay mode's weight according 
to the experimental uncertainty in its branching fraction || and the limited statistics of the 
simulation. 

The uncertainty due to the shape of the signal distribution is assessed by examining 
the variation of Xd as the data are refit with modified signal distributions. Modifications 
include variations in D*p polarization, overall tracking efficiency, solution multiplicity, and 
beam energy. The uncertainty due to the fitting distribution of combinatoric decays is 
evaluated by considering possible variations in the momentum spectrum of random tracks. 
We consider three variations, corresponding to decreasing the average momentum of soft 
tracks, increasing the average momentum of high-momentum tracks and changing the shape 
of the low-momentum spectrum as if there were 50% more D*'s than expected. 

By comparing the yield of B°B° candidate events to the yield of B° B° candidates, 
we can limit e B from this measurement ||18|| . We compare Xd in events with positively 
charged leptons to Xd with negatively charged leptons, defining x±= h w i^ / {hy V l ± + h% / l ± ) 
and Aqp = (x+ — X-)/{x+ + X-)-> where Aqp = 4-Ke(ee). Charge asymmetries in lepton 
identification cancel with this method. From studies of detection asymmetries for hadrons 
and from measurements of hadronic fake contributions, we find negligible systematic bias in 
the measurement and estimate a systematic uncertainty of 1.4% on Aqp- We determine Aqp 
= 0.017±0.070±0.014, corresponding to 9fce(e B ) = 0.004 ±0.018 ±0.003, or \&e(e B )\ < 0.034 
at the 95% confidence level. 

We are also able to provide a non-trivial limit on using Eqn.(|I|) under two assumptions. 
We assume that any possible AI^ has negligible impact upon the extraction of Arrid from 
the experimental results listed in || , and we assume that indirect CP violation is not present 
(\p/q\ = 1). We combine our measurement of Xd with the values Arrid = 0.464 ± 0.018 ps -1 
and r B o = 1.56 ± 0.04 ps |J, to find \ya\ < 0.41 at 95% confidence level. 

We have described a measurement of B° — B° mixing parameters Xd and e B , and have 
combined our result with direct measurements of Am^ to extract limits on yd. We exploit a 
method of partial reconstruction of the decay chains B — > D* ±r K T and B — > D*^p T subject 
to systematic uncertainties different from previous measurements. We note that this result 
is independent of previous CLEO mixing analyses which used leptons to tag the beauty 
quantum numbers at decay for both B mesons in the event ||. This measurement provides 
the first non-trivial limits on yd- 
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